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ABSTRACT 

The  spectrometric  oil  analysis  program  as  applied  to  naval 
aviation  was  developed  as  a  portion  of  the  overall  aviation  safety 
program  of  the  U.  S.  Navy.   The  equipment  and  techniques  have  been 
refined,  and  the  program  has  been  steadily  expanded  since  its 
inception  in  1955.   The  value  of  this  system  in  determining  densi- 
ties of  microscopic  particles  of  certain  oil-wetted  wear  metals 
in  samples  of  oil  extracted  from  aircraft  engines  has  proved  to  be 
helpful  in  predicting  incipient  engine  failure.   In  this  study 
data  relating  to  both  reciprocating  and  jet  engine  models  was  ana- 
lyzed in  an  attempt  to  determine  which  of  the  following  elements 
provided  significant  information  regarding  the  internal  condition 
of  the  engine:   aluminum,  iron,  chromium,  silver,  magnesium, 
nickel,  copper,  and  silicon. 

Multiple  and  simple  linear  regression  analyses  and  correl- 
ation techniques  were  applied  in  order  to  determine  the  mathemat- 
ical model  which  corresponded  most  closely  to  the  data  compiled. 
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I.   INTRODUCTION 

The  United  States  Navy  Spectrometric  Oil  Analysis  Program  had 
a  modest  beginning  in  19S5  at  the  Naval  Air  Rework  Facility,  Naval 
Air  Station,  Pensacola,  Florida.   From  this  beginning  the  program 
has  grown  to  the  extent  that  it  now  enjoys  a  substantial  amount  of 
credit  for  the  success  of  the  Navy's  overall  compaign  for  aviation 
safety.   The  degree  of  success  attained  was  largely  due  to  the 
ability  of  certain  types  of  spectrometers  to  detect  changes  in 
amounts  of  microscopic  particles  of  metals  suspended  in  an  oil  so- 
lution.  Experience  has  shown  that  early  identification  of  incip- 
ient engine  failure  may  be  recognized  through  proper  interpretation 
of  the  analysis  of  rates  of  particle  build-up  in  the  oil  samples. 

Current  state-of-the-art  spectrometry  can  measure  concentra- 
tions of  all  of  the  important  wear  metals  produced  in  an  oil- 
lubricated  mechanical  system.   Since  metallic  ions  emit 
characteristic  light  spectra  when  vaporized  in  an  electric  arc,  a 
unique  spectrum  for  each  metal  is  obtained.   When  the  intensity  of 
the  spectrum  is  measured  against  a  relative  standard,  the  result 
is  an  estimate  of  the  quantity  of  the  metal  present  in  the  sample 
being  analyzed. 

With  the  spectrometric  analysis  of  the  sample  completed,  the 
only  remaining  task  involves  interpretation  of  data.   The  acces- 
sion of  modern  unit  record  equipment  has  allowed  the  oil  analysis 
program  to  expand  and  to  build  a  sufficiently  large  data  base 
from  which  statistical  inferences  may  be  drawn.   This  collection 
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of  data  and  experience  led  toward  the  establishment  of  threshold 
limits  for  certain  elements  on  each  particular  engine  model. 
These  limits  have  been  delineated  for  use  as  a  guide  in  determin- 
ing whether  or  not  an  engine  is  in  need  of  overhaul. 

To  assist  engine  mechanics  in  the  performance  of  their  mainte- 
nance duties,  the  analyst  in  the  central  laboratory  at  Pensacola 
can  often  pinpoint  a  potential  trouble  spot  in  an  engine  by  noting 
which  particular  metal  has  exceeded  its  threshold.   For  instance, 
aluminum  is  often  used  in  reciprocating  engines  for  parts  such  as 
pistons,  cylinder  heads,  and  oil  pump  impellers.   Copper  is 
almost  always  used  for  bushings  and  intake  valve  guides,  and  nickel 
is  used  for  exhaust  valve  guides.   No  single  standard  list  can  be- 
compiled  which  would  be  appropriate  for  all  engines.   However, 
each  model  has  a  list  of  standardized  parts.   This  standardization 
within  models  is  what  allows  the  analyst  at  the  central  laboratory 
to  detect  not  only  probable,  impending  trouble  but  also  to  predict 
the  engine  part  which  will  most  likely  fail. 
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II.   CONCEPT  OF  THE  STUDY 

This  study  will  be  concerned  primarily  with  analysis  of  data 
collected  during  the  period  1  July  1967  -  30  September  1967  at  the 
Naval  Air  Station,  Pensacola,  Florida.   Attention  will  be  focused 
upon  a  reciprocating  engine  model  produced  by  the  Pratt  and  Whitney 
Aircraft  Division,  United  Aircraft  Corporation.   The  specific  model 
is  the  R2000-6.   The  data  for  three  different  models  of  jet  engines 
was  also  examined.   The  jet  engine  analysis  covered  models 
JS2-P-6,  J60-P-6,  and  J33-A-24.   Since  the  results  relating  to  jet 
engine  data  differed  from  those  relating  to  the  reciprocating  model 
a  discussion  of  general  results  for  the  jet  models  is  presented 
separately  in  Appendix  H.   The  procedure  for  spectrometric  analysis 
of  oil  samples  is  identical  for  both  jet  and  reciprocating  engines. 
Whenever  a  sample  is  analyzed  by  the  spectrometer,  the  data  is 
automatically  recorded.   Each  sample  provides  a  data  point  consist- 
ing of  a  number  of  components .   Fifteen  of  the  components  which  are 
of  interest  to  this  study  are 

(1)  Engine  model  number 

(2)  Engine  serial  number 

(3)  Date  oil  sample  removed  from  engine 

(MO  Hours  since  engine  overhaul  (engine  operating  hours) 

(5)  Hours  since  engine  oil  change  (engine  operating  hours) 

(6)  Density  of  aluminum  (parts  per  million) 

(7)  Density  of  iron  (parts  per  million) 

(8)  Density  of  chromium  (parts  per  million) 
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(9)  Density  of  silver  (parts  per  million) 

(10)  Density  of  copper  (parts  per  million) 

(11)  Density  of  tin  (parts  per  million) 

(12)  Density  of  magnesium  (parts  per  million) 

(13)  Density  of  lead  (parts  per  million) 

(14)  Density  of  nickel  (parts  per  million) 

(15)  Density  of  silicon  (parts  per  million) 

Throughout  this  study  we  shall  refer  to  several  of  the  above- 
mentioned  elements  less  frequently  than  to  others.   For  instance, 
silicon  is  the  only  non-metallic  element  in  the  data.   It  is  used 
primarily  as  an  indicator  of  the  amount  of  sand  that  a  reciprocat- 
ing engine  has  ingested  from  the  atmosphere.   Analysis  for  silicon 
contamination  in  jet  engines  has  become  meaningless  due  to  the 
increased  use  of  silicone  additives  in  jet  engine  oils. 

Reference  to  a  specific  serial  number  of  a  specific  model 
engine  will  be  by  the  digit  numbers  (1) ,  (2) , . . . ,  (26)  rather  than 
by  the  full  six  digit  serial  numbers.   This  simplifies  the  nota- 
tion throughout  the  study.   Appendix  A  may  be  referred  to  for 
determining  which  serial  number  is  related  to  which  digit  number. 

Practically  speaking,  the  variables  which  can  be  measured 
without  error  are  the  hours  since  engine  overhaul  and  the  hours 
since  engine  oil  change.   These  two  variables  have  been  designated 
the  independent  variables,  and  the  densities  of  the  ten  elements 
have  been  designated  as  the  dependent  variables.   The  purpose  of 
this  study  is  the  investigation  of  relationships  which  may  exist 
between  items  (4)  and  (5)  of  the  list  of  components,  and  the 
last  ten  items. 
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Since  little  information  of  an  analytical  nature  exists,  the 
observed  data  was  plotted  on  a  scatter  diagram  to  aid  in  determin- 
ing what,  if  any,  functional  relationships  exist.   With  acceptance 
of  the  assumption  of  normality,  regression  analysis  provided  the 
natural  technique  for  analytically  uncovering  linear  trends 
between  the  variables,  and  for  expressing  any  existing  functional 
relationships  in  mathematical  form.   By  this  method  the  best  rela- 
tion among  the  variables  was  determined.   Further  analysis  using 
correlation  techniques  showed  the  strength  with  which  the  variables 
were  associated.   When  preliminary  investigations  into  the  data 
showed  a  rather  high  correlation  between  the  two  independent  vari- 
ables, a  stepwise  regression  procedure  was  conducted  in  order  to 
determine  the  more  desirable  course  of  action:   continuation  of 
the  analysis  using  multiple  linear  regression  techniques,  or 
simplifying  the  analysis  by  employing  simple  linear  regression. 
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III.   PROCEDURE 

Computations  and  computer  plots  for  this  study  were  conducted 
on  the  IBM  360  computer  at  the  Naval  Postgraduate  School, 
Monterey.   Initial  scatter  diagrams  were  not  encouraging.   The 
plot  of  each  element  versus  each  of  the  two  independent  variables 
was  reproduced  in  Appendix  B.   These  plots  depicted  the  entire  set 
of  data  points  for  all  26  engines  for  model  R2000-6.   Although  the 
data  did  not  appear  to  have  strong  linear  properties,  the  plots 
were  interesting  enough  to  invite  a  more  detailed  investigation, 
particularly  the  trends  exhibited  by  aluminum,  iron,  silver,  and 
copper. 

A  sampling  of  five  of  the  engines  was  selected  for  individual 
plots  and  the  data  thus  obtained  graphed  separately.  The  results 
are  shown  in  Appendix  C.  These  results  proved  to  be  far  more 
encouraging  since  they  reinforced  the  initial  impression  that  the 
four  elements;  aluminum,  iron,  copper,  and  silver  showed  definite 
linear  trends. 

Multiple  linear  regression  techniques  were  applied  to  the  raw 
data.   The  predicted  value,  (yT) ,  was  expressed  by  the  equation 

y1  =  a  +  b^  xi  +  b2  X2 
Here  x,  is  hours  since  engine  overhaul;  and  x^  is  hours  since 
engine  oil  change.   These  are  the  variables  affecting  y  which  is 
the  density  in  parts  per  million  of  the  element  being  predicted. 
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Estimates  for  the  partial  regression  coefficients,  b^  and  02,  and 
the  constant  term,  a,  were  computed  using  the  procedure  outlined 
in  Ostle. 

While  analyzing  the  26  engines  of  model  R2000-6,  it  was  consid- 
ered desirable  to  determine  if  all  26  engines  could  be  treated  as 
if  they  came  from  the  same  population.   In  addition  to  the  normal- 
ity assumption,  if  the  further  assumption  of  equal  variances  for 
all  26  engines  proved  acceptable,  the  data  could  possibly  be  pooled 
and  the  engines  could  be  handled  as  one.   However,  due  to  the 
suspect  nature  of  the  initial  scatter  diagrams,  the  assumption  of 
homogeneity  of  variances  was  not  made.   Therefore,  the  test  selec- 
ted  to  handle  this  problem  was  the  one  developed  by  Bartlett.    In 
order  to  test  for  homogeneity  of  variances  the  null  hypothesis  was 

h    cr 2  =  G~2  =         =CT2=  rr2 

nQ  .   vj1     \-i2     ...     vj26     <j 
Bartlett' s  statistic  has  been  shown  to  be  distributed  under  HQ 
approximately  as  chi  square  with  K-l  degrees  of  freedom. 

The  test  was  conducted  at  the  OC  =  .05  level  for  the  elements 
aluminum,  iron,  chromium,  silver,  magnesium,  nickel,  copper,  and 
silicon.   With  the  exception  of  chromium,  the  null  hypothesis  was 
rejected.   Not  only  had  the  equal  variances  test  failed,  but  it 
was  also  found  that  the  multiple  linear  regression  had  indicated 
that  the  partial  regression  coefficient  b^  would  often  have  a 


-'-Ostle,  B.s  Statistics  In  Research,  2nd  Ed.,  p.  177,  Iowa 
State,  1963. 

2 
Brownlee,  K.  A.,  Statistical  Theory  and  Methodology  In 

Science  And  Engineering,  p.  225,  Wiley,  1960. 
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positive  value  while  b2  would  have  a  negative  value  and  vice  versa, 
Since  the  evidence  was  overwhelmingly  against  homogeneity  of  vari- 
ances for  the  multiple  linear  regression  situation,  it  was  deemed 
advisable  to  conduct  a  stepwise  regression  analysis  to  determine 
which  of  the  two  independent  variables  was  the  better  predictor. 
A  further  benefit  from  this  analysis  was  that  an  objective  evalu- 
tion  could  also  be  made  to  determine  the  relative  merit  of  retain- 
ing the  other  variable  in  the  equation. 

The  stepwise  procedure  used  was  a  linear  regression  program 
designated  BIMED02R.   This  program  is  one  of  a  series  of  library 
programs  designed  by  the  Health  Sciences  Computing  Facility, 
University  of  California,  Los  Angeles.   A  sequence  of  multiple 
linear  regression  equations  is  computed  in  a  stepwise  manner.   At 
each  step  one  variable  is  added  to  the  regression  equation.   The 
variable  added  is  the  one  which  makes  the  greatest  reduction  in 
the  error  sum  of  squares.   This  is  the  same  as  saying  that  it  is 
the  variable  which  has  the  highest  partial  correlation  with  the 
dependent  variable  partialed  on  the  variables  which  have  already 
been  added  to  the  equation.   Since  this  analysis  had  only  two  inde- 
pendent variables,  the  first  variable  which  was  entered  into  the 
equation  "explained  more  of  the  variance"  than  did  the  second 
variable . 

Stepwise  regressions  were  conducted  on  each  of  the  elements 
for  which  the  Naval  Air  Rework  Facility  had  previously  established 
threshold  limits.   For  engine  model  R2000-6  limits  have  been 
established  for  aluminum,  iron,  chromium,  silver,  copper,  nickel, 
and  silicon.   Results  of  the  stepwise  regression  for  the  four 
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elements  aluminum,  iron,  copper,  and  silver  showed  that  hours  since 
engine  oil  change  appeared  as  the  first  variable  entered  in  the 
regression  equation  far  more  often  than  did  hours  since  engine 
overhaul.   The  results  obtained  from  the  remaining  elements  were 
inconclusive  as  to  which  of  the  two  independent  variables  was  more 
consistent  in  explaining  variance.   For  these  elements  hours  since 
engine  overhaul  appeared  as  the  first  variable  entered  in  the  step- 
wise procedure  about  as  often  as  did  hours  since  engine  oil  change . 
Appendix  D  contains  the  specific  outcome  for  each  of  the  26 
engines . 

With  the  stepwise  regression  completed,  a  test  of  the  hypothe- 
sis that  the  slope  of  the  regression  line  was  zero  was  conducted. 

The  first  test  had  as  the  null  hypothesis: 

H01  :   9li  =  °     ;  i  =  1,  2,...,  26 
To  accomplish  this,  the  following  ratio  was  formed  from  the  Anal- 
ysis of  Variance  Table: 

mean  square  due  to  bjla 
residual  mean  square 

This  ratio   follows  an  F-  distribution  with   M,  =  1  and 
\J   p  =  ni  _  2  degrees  of  freedom  where  n-^  is  the  number  of  obser- 
vations taken  on  the  i-th  engine. 

The  next  test  had  as  the  hypothesis: 

H02  :   ^2i  =  °     i1'1-  2"-"  26 
To  accomplish  this  test,  the  following  ratio  was  formed  from  the 
Analysis  of  Variance  Table: 

ian  square  due  to  b2|a,  b, 
residual  mean  square 
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This  ratio  follows  the  F-  distribution  with  P -,=  1  and  M?  =  n±~^ 
degrees  of  freedom.  The  results  of  testing  the  foregoing  hypothe- 
sis at  the  0C=  .05  level  are  contained  in  Appendix  E. 

The  ratio  for  testing  Hm  assesses  the  significance  of  the 
reduction  in  the  residual  sum  of  squares  gained  by  entering  the 
first  variable  in  the  equation.   The  notation  bjla  is  used  to  indi- 
cate that  the  sums  of  squares  associated  with  the  coefficients  are 
obtained  sequentially,   Likewise,  the  ratio  for  testing  Hq2  is 
used  to  determine  if  any  significance  should  be  attached  to  the 
additional  reduction  in  the  residual  sum  of  squares  achieved  by 
entering  b2  into  the  equation  after  the  sums  of  squares  for  a  and 
b-L  have  been  computed.   Thus,  the  notation  b2  a,b^  was  used. 

As  part  of  the  stepwise  procedure  the  estimated  multiple  cor- 
relation coefficients  were  computed  at  the  completion  of  each 
step.   The  addition  of  the  second  variable  into  the  multiple  lin- 
ear regression  equation  caused  an  increase  in  the  correlation 
coefficient  of  0.1000  or  more  in  only  the  number  of  engines  shown 
in  Table  I. 

The  relatively  small  increase  in  correlation  is  not  considered 
significant  in  view  of  the  additional  increase  in  the  computational 
effort  required  to  effect  the  higher  correlation. 
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Element 

Number 

of  Engines 

Aluminum 

4 

Iron 

2 

Copper 

4 

Silver 

5 

Chromium 

3 

Magnesium 

9 

Nickel 

7 

Silicon 

6 

TABLE  I 


Increases  of  0.10  or  More  in  Multiple  Correlation 
Coefficient  Caused  by  Addition  of  the  Second 
Independent  Variable  to  the  Equation 
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IV.   RESULTS 

A  significant  result  obtained  from  the  analysis  of  the  step- 
wise regression  showed  that  hours  since  engine  oil  change  was  the 
better  predictor  of  the  two  independent  variables.   Furthermore 
when  the  results  of  the  tests  for  zero-slope  were  scanned,  it  was 
noted  that  the  hypothesis 

H02  :   ^2i  =  °    ;  i  =  1,  2,...,  26 
was  accepted  more  often  than  the  hypothesis 

H01  :   9li  "  °    ;  i  =  1,  2,...,  26 
This  acceptance  of  hypotheses  held  true  for  each  of  the  metallic 
elements  in  the  analysis.   For  the  element  silicon  the  number  of 
engines  for  which  Hq2  was  accepted  proved  to  be  equal  to  the  num- 
ber for  which  Hqq_  was  accepted.   This  result,  when  considered  with 
the  multiple  correlation  coefficient  results,  indicated  that  not 
only  did  the  addition  of  the  second  variable  to  the  regression  add 
little  to  the  correlation  coefficient  but,  in  some  cases,  its 
addition  actually  tended  to  nullify  some  of  the  effect  of  positive 
slope  provided  by  the  first  variable.   That  is,  it  was  found  that 
instances  occurred  wherein  one  independent  variable  was  positively 
correlated  with  the  dependent  variable,  and  the  other  independent 
variable  was  negatively  correlated  with  the  dependent  variable. 
With  this  evidence  established  the  decision  was  made  to  drop  the 
variable  hours  since  engine  overhaul  from  the  analysis  and  to 
proceed  using  only  the  remaining  independent  variable  hours  since 
engine  oil  change. 
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The  computer  program  in  Appendix  G  was  written  to  perform  a 
simple  linear  regression  on  each  of  the  26  engines  for  each  of  the 
following  elements:   aluminum,  iron,  chromium,  silver,  magnesium, 
nickel,  copper,  and  silicon.   The  possibility  of  constructing  one 
single  regression  equation  which  could  be  used  for  all  26  engines 
had  already  been  ruled  out  by  Bartlett's  Test  for  homogeneity  of 
variances.   Therefore,  the  conclusion  that  each  engine  should  be 
treated  individually  was  firmly  established. 

A  final  task  which  the  program  in  Appendix  G  was  designed  to 
perform  was  that  of  establishing  a  prediction  interval  of  the  inde- 
pendent variable  x.   Such  an  interval  was  to  be  predicted  for  each 
engine  and  for  each  element. 

Now,  since  the  regression  was  reduced  to  a  simple  linear  form, 
the  expected  value  of  y  which  corresponds  to  a  given  x  =  X  is 
calculated  as 

yT  =  a  +  bX 

Using  this  predictor  for  y  we  can  employ  the  formula  for  a 

3 
100  (1-0C)%  prediction  interval  for  y  which  is  presented  in  Ostle. 

Appendix  F  displays  the  results  obtained  from  a  hypothetical 

value  of  x  =  150  for  a  95  per  cent  prediction  interval  around  y. 

Since  y  represented  the  parts  per  million  density  for  the  element 

under  consideration,  the  low  end  of  the  interval  was  set  equal  to 

zero  at  any  time  that  the  predicted  low  end  dropped  below  zero. 

A  negative  number  of  parts  per  million  for  an  element  would  be 


30stle,  B.,  Statistics  In  Research,  2nd  Ed.,  p.  170,  Iowa 
State,  1963. 
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meaningless.   Also,  since  the  spectrometer  in  use  at  the  Naval  Air 
Station,  Pensacola,  measures  densities  of  the  elements  in  integer 
values,  the  low  end  values  were  truncated  to  integer  values  and  the 
high  end  values  were  rounded  up  to  the  next  highest  integer  if  they 
were  composed  of  any  decimal  fraction.   An  illustration  is  shown  in 
Drawing  I . 
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V.   CONCLUSIONS 

The  policy  pertaining  to  engine  model  R2000-6  requires  col- 
lection of  data  for  both  of  the  independent  variables  hours  since 
engine  oil  change  as  well  as  hours  since  engine  overhaul.   In  ad- 
dition to  this,  data  is  also  collected  for  aluminum,  iron,  chromi- 
um, silver,  copper,  nickel,  and  silicon.   It  appears  that 
virtually  all  of  the  information  relating  to  this  model  engine  is 
contained  in  the  data  for  aluminum,  iron,  copper,  silver,  and 
hours  since  engine  oil  change.   These  five  items  provide  a  notice- 
ably linear  trend  whereas  chromium,  nickel,  silicon,  and  hours 
since  engine  overhaul  appear  to  possess  poorer  predictive  powers. 

Each  engine  has  its  own  characteristic  build-up  rate  for  each 
element;  therefore  data  collection,  storage,  and  manipulation 
must,  of  necessity,  be  accomplished  separately  for  each  particular 
serial  number  of  each  model  engine.   When  a  new  engine  serial 
number  is  introduced  into  the  system,  there  can  be  very  little 
information  of  a  predictive  nature  which  can  be  fed  back  to  the 
operating  unit  until  a  sufficient  number  of  data  points  has  been 
recorded.   This  study  utilized  serially  numbered  engines  only  if 
there  were  eight  or  more  data  points  with  which  to  work.   Once  an 
engine  has  provided  enough  data  points  to  characterize  its  build- 
up rate  of  parts  per  million  for  each  element,  any  new  reading  can 
bo  examined  against  its  prediction  interval  and  appropriate  action 
taken.   Normally,  appropriate  action  could  be 

(1)  If  the  actual  y  value  is  within  the  prediction  interval  for 
all  elements,  the  point  is  added  to  the  data  bank  and  the 
prediction  mechanism  uses  this  data  point  along  with  all 
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previous  points  to  update  the  regression  equation.   Obviously, 
the  more  data  points  collected  on  a  particular  engine,  the 
better  should  be  the  prediction. 

(2)    If  the  actual  y  value  is  above  the  predicted  high  end  of  the 
interval  for  any  one  or  more  of  the  elements  associated  with 
a  particular  model  engine,  then  the  operating  unit  is  noti- 
fied with  an  appropriate  recommendation.   The  data  point  is 
then  added  to  the  data  bank  and  treated  as  in  (1)  above. 

A  further  study  should  be  conducted  using  data  collected  over 

a  longer  period  of  time  such  as  a  three-year  period.   The  results 

obtained  from  such  a  study  should  be  much  more  conclusive  than 

these  findings  which  have  been  based  on  data  collected  over  a 

three-month  period. 
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APPENDIX  A 

SERIAL  NUMBERS  OF  THE  R20Q0-6  MODEL  AIRCRAFT  ENGINES 

Engine  Number  Serial  Number 

1  101196 

2  103105 

3  104166 

4  105005 

5  105246 

6  106505 

7  106666 

8  106681 

9  106690 

10  106707 

11  106797 

12  106802 

13  107486 

14  107935 

15  108003 

16  108372 

17  108458 

18  108461 

19  108553 

20  108616 

21  108919 

22  700002 

23  70013u 

24  700163 

25  702269 

26  702489 
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APPENDIX  C 


Plots  of  Hours  Since  Engine  Oil  Change 
Versus  Densities  of  Eight  Elements  for  Five  Engines 
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APPENDIX  D 
STEPWISE  REGRESSION  RESULTS 


First  Indepe 

;ndent 

Vari. 

able  Ente 

red 

in  Stepwise 

Regression 

Engine  No. 

Mg 

Al 

Fe 

Cu 

Ag 

Cr 

Ni 

Si 

1 

2 

2 

2 

1 

2 

2 

2 

2 

2 

1 

2 

2 

2 

1 

2 

2 

1 

3 

2 

2 

2 

2 

2 

1 

1 

2 

4 

2 

2 

2 

2 

2 

2 

1 

2 

5 

2 

2 

2 

2 

2 

1 

2 

1 

6 

1 

1 

1 

2 

2 

2 

1 

2 

7 

1 

2 

2 

2 

2 

2 

1 

1 

8 

2 

2 

1 

2 

2 

1 

1 

2 

9 

2 

2 

2 

2 

2 

1 

1 

1 

10 

1 

1 

2 

2 

2 

2 

2 

2 

11 

2 

2 

1 

2 

1 

2 

1 

2 

12 

1 

2 

2 

2 

2 

1 

1 

1 

13 

2 

2 

2 

2 

2 

1 

2 

1 

14 

1 

1 

1 

1 

1 

2 

1 

2 

15 

1 

1 

1 

1 

1 

1 

2 

2 

16 

2 

2 

2 

2 

2 

2 

2 

1 

17 

2 

1 

2 

2 

2 

1 

2 

2 

18 

2 

2 

2 

2 

2 

1 

2 

2 

19 

2 

2 

2 

2 

2 

1 

1 

1 

20 

2 

1 

2 

2 

2 

2 

2 

2 

21 

1 

1 

2 

2 

2 

1 

2 

2 

22 

1 

2 

2 

2 

2 

2 

1 

2 

23 

1 

1 

2 

2 

2 

1 

1 

2 

24 

2 

2 

2 

2 

2 

1 

1 

2 

25 

1 

1 

1 

2 

1 

1 

1 

1 

26 

2 

2 

2 

2 

1 

1 

2 

2 

Total  1 

11 

9 

6 

3 

6 

15 

14 

9 

Total  2 

15 

17 

20 

23 

20 

11 

12 

17 

Key:   1   Signifies  that  the  independent  variable,  hours  since 
engine  overhaul,  was  the  first  variable  entered  in 
the  stepwise  regression. 

2   Signifies  that  hours  since  engine  oil  change  was  the 
first  variable  entered  in  the  stepwise  regression. 
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APPENDIX  F 


DATA  FOP  ALUMINUM 


PREDICTION  INTERVAL  FOR  AN  INDIVIDUAL  VALUEOF  Y  FDR  A 
HYPOTHETICAL  VALUE  OF  THE  INDEPENDENT  VARIABLE* X  =  150 


PREDICTED  PREDICTED 

ENGINE  NO.  LOW  END  HIGH  END 

1  3.00000  27.00000 

2  0.0  10.00000 

3  R. 00000  27.G00C0 
^  0.0  19.00000 

5  4.00000  32.00000 

6  0.0  39.00000 

7  1.00000  24.00000 

8  3.00000  21.OG0CO 

9  2.00000  27.00000 
1C  2.00C00  17.00000 

11  4.03000  25.00000 

12  3.00000  53.00000 

13  11.00000  29.00000 

14  8.00000  30.00000 

15  0.0  33.00000 

16  0.0  22.00000 

17  14.00000  23.00000 

18  6.00000  20.00000 

19  3.00000  20.00000 

20  3.00000  18.00000 

21  2.00000  44.00000 

22  0.0  19.00000 

23  5.00000  28.00000 

24  6.00000  13.00000 

25  2.00000  32.00000 

26  4.00000  11.00000 
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DATA  FOR  IRON 


PREDICTTCN  INTERVAL  FOR  AN  INDIVIDUAL  VALUEOF  Y  FOR  A 
HYPOTHETICAL  VALUE  OF  THE  INDEPENDENT  VARIABLE, X  =  15C 


ENGINE  NO. 


PREDICTED 
LOW  END 


PREDICTED 
HIGH  END 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 


49.00000 
25.00000 
47.00000 
45.00000 
59.00000 
19.00000 
38.00000 
91.00000 
40.00000 
38.00000 
47.00000 
62.00000 
56.00000 
32.00000 
30.00000 
17.00000 
79.00000 
38.00000 
44.00000 
44.00000 
42.00000 
44.00000 
35.00000 
30.00000 
26.00000 
29.00000 


95. 

54. 

82. 

73. 
114. 
123. 
151. 
166. 
105. 

72. 

87. 
124. 
109. 
119. 
118. 

77. 
102. 

76. 

89. 

85. 
147. 

80. 
112. 

66. 
136. 

53. 


00000 
OOOCO 
00000 
00000 
00000 
00000 

coooo 

00000 
00000 
00000 
OOOCO 
00000 
00000 
00000 
00000 
00000 
00000 
00000 
00000 
00000 
00000 
00000 
OOOCO 
00000 
00000 
00000 
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DATA    FOR    COPPER 


PPFDICTICN    INTERVAL    FOR    AN    INDIVIDUAL    VALUEOF    Y    FOR    A 
HYPOTHETICAL    VALUE    OF    THE    INDEPENDENT    VARIA8LE,X    =     15C 


ENGINE  NO, 


PREDICTED 
LOW  END 


PREDICTED 
HIGH  END 


1 
2 
3 

u 

5 

6 

7 

8 

Q 

10 

11 

12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 


6.0J0Q0 
4.00000 

10.00000 
7.00000 

10.00000 
0.0 
4.00000 

10.00000 
4.00000 
4.00000 
5.00000 
9.00000 
9.00000 
0.0 

5.00000 
3.00000 

18.00000 
8.00000 
0.0 

6.00000 
7.00000 
1.00000 
2.00000 
8.00000 
2.00000 
5.00000 


30.000C0 
13.0000C 
18.00000 
20.00000 
27.000C0 
22.00000 
26.00000 
25.00000 
26.00000 
14.00000 
19.00000 
19.00000 
27.00000 
30.00000 
20.00000 
21.00000 
27.00000 
18.00000 
23.00000 
20.00000 
23.000C0 
19.00000 
13.00000 
22.00000 
32.00000 
17.00000 


i*8 


OAT*     i=TR    flLVER 


PPEOICTTCN    INTERVAL    FPP-     AN    INOIVIDUAL    VALMFOF    Y    FOP     A 
HYPGTHFTICAL    VALUE    "F    THE     INDEPcNOENT    VARIABLE, X    =     L5C 


PPECICTEC  PREOICTFn 

FNGINE     NO.  LOW    ENO  HIGH    END 

1  0.0  12.000C0 

2  0.0  5.000CO 

3  A. 00000  14.0COCO 
*  1.00000  12.0J0OO 

5  3.0JO00  16.00000 

6  0.0  18.OO0C0 

7  0.0  18.000CO 

8  0.0  11.00000 

9  0.0  16.C00C0 
1C  l.OOCOO  9.00000 

11  1.00000  8.000CO 

12  0.0  10.00000 

13  4.0000  0  12.000CJ 

14  I. 00000  16.C00CO 

15  0.0  13.00000 

16  0.0  ll.OOOCO 

17  3.00000  14.00000 

18  4.00000  14.00000 

19  0.0  13.00000 

20  1.00000  13.00000 

21  1.00000  ll.OOOCO 

22  0.0  9.00000 

23  0.0  7.00000 

24  0.0  1C.C00C0 

25  0.0  18.00000 

26  0.0  8.00000 


«+9 


DATA     FOP    CHROMIUM 


PPFHTC'ICN    INTERVAL    FOP    AN    INDIVIDUAL    VALUEOF    Y    FOR    A 
HYPOTHETICAL    VALUF    OF    THE    INDEPENDENT    VARIABLE, X    =     15C 


F  NIG  INF  NO. 


1 
2 
3 
A 
* 
6 

7 

P 
Q 

10 
11 

12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 


PREriCTED 

LOW  END 

0 

.0 

0 

.0 

0 

.0 

0 

.0 

0 

.0 

1 

.00000 

0 

.0 

0, 

.0 

0 

.0 

0 

.0 

0 

.0 

0 

.0 

1 

.00000 

0 

.0 

0 

.0 

0 

.0 

1 

.00000 

0 

.0 

0 

.0 

0 

.0 

2 

.00000 

0 

.0 

0 

.0 

0 

.0 

0 

.0 

0 

.0 

PREDICTED 
HIGH  END 


19.00000 
12.000C0 
16.00000 
16.00000 
18.000C0 
21.00000 
24.00000 
2  3.OO0OC 
20.00000 
18.00000 
39.00000 
27.00000 
21.00000 
25.00000 
25.00000 
17.00000 
25.00000 
13.00000 
16.00000 
14.00000 
30.00000 
13.00000 
23.00000 
16.00000 
18.00000 
15.000C0 
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DATA    ?fip    magncSTUm 


PPFHICTICN    INTERVAL    FOR    AN    INDIVIDUAL    VALUEOF    Y    Fne     A 
HYPOTHETICAL    VALUF     HF    ^HF     INDFPcNnFNT    V*RURLE,X    =     1  «5 r 


F  NO  INF     NO, 


1 
2 
3 
i. 

5 

6 
7 
B 

ir 
11 
12 
13 
1* 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 


PPE 

niCTFC 

LOW  FNO 

0 

.0 

0 

.0 

2 

.00000 

0 

•0 

0 

.0 

0 

.0 

0 

.0 

0 

.0 

0 

.0 

0 

.0 

0 

.0 

0 

.0 

3 

.00000 

0 

.0 

0 

.0 

0 

.0 

1 

.00000 

0 

.0 

0 

.0 

0, 

,c 

0, 

.0 

0 

.0 

0 

.0 

0. 

.0 

0 

.0 

0 

.0 

PREOICTFP 

HIGH  FMO 

4, 

,00000 

2, 

.00000 

6, 

,00000 

3, 

,00000 

6, 

,00000 

7, 

,00000 

4, 

.OOOGC 

6, 

(    -   V-   J  »,   ./ 

6, 

,00000 

2, 

.oocco 

4, 

,0000  0 

4, 

,00000 

11, 

,00000 

6, 

,00  000 

6, 

,00000 

3, 

,00000 

5, 

.CCGOU 

4, 

.COOCO 

3, 

,00000 

3, 

,00000 

5, 

,00000 

4, 

,0000  0 

6. 

,00000 

3, 

,oooco 

7. 

,00000 

3, 

,00000 
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DATA    FOR    NICKEL 


PPEDICTICN    INTERVAL    FOR    AN    INDIVIDUAL    VALUEOF    Y    FOR     A 
HYPOTHETICAL    VALUE    OF    THE    INDEPENDENT    VARIABLE, X    =     ISC 


PREDICTED  PREDICTED 

rNGINF    NO.  LOW    END  HIOH    END 

1  0.0  7.000'JC 

2  0.0  3.00000 

3  0.0  4.0C0C0 

4  0.0  A. 00000 

5  0.0  5.00000 

6  0.0  5.00GC0 

7  0.0  6.00000 

8  0.0  5.00000 

9  0.0  5.00C00 

10  0.0  6.000CC 

11  0.0  4.000C0 

12  0.0  3.000C0 

13  0.0  5.00000 

14  0.0  5.000CO 

15  J.O  5.00000 

16  0.0  4.00000 

17  0.0  3.000CO 

18  0.0  3.00000 
1Q  0.0  3.00000 

20  0.0  5.00000 

21  0.0  6.00000 

22  0.0  4.00000 

23  0.0  5.C00C0 
2^  0.0  4.00000 

25  0.0  8.00000 

26  0.0  3.OO0C0 
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DATA  *0P  MLICCN 


PREDICTTCN  INTFRV^L  FOP  AN  I  NO  I V  !DU  AL  vaLHFOF  Y  FOR  A 
HYPCTHFTICAL  VALUE  OF  "-"HE  INDFPcNOFNT  VARIABLE, X  =  15C 


FNGINE  NO 

1 

2 

3 

4 

5 

6 

7 

8 

o 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

PREDICTEC 
LOW  END 


O.C 

0.0 

1.00000 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

l.COCOO 

1.00003 

0.0 

1.00000 

0.0 

0.0 

0.0 

1.00000 

0.0 

0.0 

0.0 

0.0 

0.0 


PRECICTFD 
HIGH  EMO 


10.00  000 

7.OO0C0 

8.000CO 

6.00000 

31.00000 

10.00000 

8.00000 

9.C00C3 

11.00000 

22.C00C0 

8.00000 

18.00000 

13.0 COCO 

7.00000 

8.000CC 

7.00000 

6.00000 

7.00  000 

1C. 00000 

7.00000 

8.00000 

7.000C0 

9.00000 

8.CO0CO 

9.00000 

7.00000 


53 


APPENDIX  G 


Computer  Program 


C 

r 

c 

c 

r 

r 


THI«  P 
AND  PE 


P^AL  * 

RFAL  * 
17,15)  , 
2)  , NIC  2 

INTEGE 

DATA 

1*0.0/, 
2105246 
3  '10670 
4,'1C80 

6     ',  '70 


RCGRAM    TAKES    THE    INDEPENDENT    VARIAEL 
GPESSCS    IT    AGAINST    THE    DEPENDENT    VAR 


no    i    i 
READ! 5 
lc14 

7     PRPMA^ 

no    I* 
IF    (     D 
14    COMTIN 

GO  TO 
31  K  L  1=  I 
LA  =  I 
DATE(L 
HOV(L , 
HOC  f  L  , 
*L(L,L 
FE(L,l 
CR(L,L 
AG(L,L 
CU(t,L 
MGCL,L 
NICL,L 
SICL,L 
1  C.PNTIN 
WPITC 

2CC  FHPMAT 
CALL  R 
WR ! TF  ( 

201  «=QPMAT 
CALL  R 
WPI TEC 

2C2    CHPMAT 

CALL    R 
WRITF( 

2C3  FORMAT 
CALL  R 
WPITEC 

2C6  FORMAT 
CALL  R 
WPITE  C 

20  5  POP  MAT 
CALL  R 
WITFC 

2C<  FOPMftj 
CALL  R 
WRITFC 

?07    POPMAT 
CALL    R 
STOP 
END 


«  AMQD/'P 
*  HPVC27, 
A0C27,15) 
7, 15), SIC 
R  *  2  K 
HO V, HOC, A 
SN/»10119 
• ,»1C65 
7  »,'106 
03  ','IC 
616  »,»1 
2269  »,» 
X=l,241 
,7)    A,D,F 


20C06 
15) .HOC 
,CUC27, 
27,15) , 
27) /2^ 

L,FE,CR 
6   »,  «1 

C"?   •,  • 

797   •  , 
8372   • 
08919 
702489 


,/,SN(27)  ,0,A,STPRFC 
C27,15)  ,ALC27,15)  ,FE 
15) ,TNC27  ,15)  ,MGC27, 
DATEC27,15) 
♦0/ 

,AG,CU,TN,MG,TI,NI,S 
•104166  •, 
,'106681 


03105 
106666 
•106802 
,•108458 
•  ,»700002 


•,•105189 


•  ,  •  107486 
•,•108461 
•  ,'70C 
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E,OIL    CHANGE 
IABLC,     PPM 


1C  )/10*0.G/ 
C27,15),CP(2 
1*)  ,TI(27, 15 


I. DATE  Z5265 
•10  500  5  •,• 
, • 106690  • , 
•,•107935  • 
•,  '108553 
•  , •700163 


C1X,A 
L  =  l,2 
.EO. 
UF 
1 

CI  )+l 
(L) 
,LA  ) 
LA)     = 
LA)    = 
A)    = 
A)    = 
A)     = 
A  )    = 
A)     = 
A  )    = 
A)     = 
A)=El 
L'E 

6,200 
CIS 
EGPFS 
6,201 
CMS 
EGPFS 
6,202 
C  »1». 
EGPES 
6,203 
CIS 
EGPES 
6,206 

(  'I'l 

EGRES 

6,205 
CI', 
EGPES 
S204 
C  »1  •, 
FGRES 
6,207 
C  «1  •. 
FGRES 


2,E3,E4,E5,E6,E7,E8,E9,E1C,E11,E12,E13, 

10F3.C  ) 


8,1X,A8,14X,A4,2X,F^.C,F3.0,5X, 
6 
SNCD)    GO    TO   31 


=    E2 

E3 
E4 

E5 

E6 

E7 

E6 

E9 

Ell 

E13 

4 

) 


///57X, 

(HOC  ,AL 

) 

///57X, 

CHOC ,FE 

) 

///57X, 

CHOCCP 

) 

///57X 

C 

) 


HOC  ,AG 


///57X, 

( HOC , MG 

) 

///57X, 

CHOC ,NI 

) 

///*7X, 

(HOC  ,CU 

) 

///*7X, 

CHOC f SI 


DATA    FOR    ALUMINUM' ) 
Itl) 

DATA    FOR     IRON') 
1.2) 

DATA    FOR    CHPCMIUM« ) 
I  .3) 

DATA    FOR    SILVER' ) 
1,4) 

DATA    FOR    MAGNESIUM') 
I  ,5) 

DATA  FOR  NICKEL') 
I  ,6) 

DATA  FOR  COPPER') 
I  .7) 

DATA  FOR  SILICON* ) 
I  ,8) 
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SIJRcntjTTME    PEGRFS(  X,Y,I  tNI)VEL) 

THIS    SUBPROGRAM    USFS     X    AS    T"F       TNDEPENDFNT    VARIABLE 
AND    Y    TS    THE    DFPFN^ENT    V^RI^8LC.    THP    I     IN   THE 
APHIMPNT     IS    THF     NUM8ER     OF     C8^c  R  VAT  I ON*     PN    F  4CH 
SEP  TALLY    NUMBERED    "200.06    fNGINc 

DIMENSION  X( 27,17)  ,  Y( 27  ,17) , I< 27 ) , ALFA ( 27 ) ,XMEAN ( 27 ) 

l,CAPY(27t 17) ,  SSQP(27)  ,Y WEAN (27)  ,PETM27 ) 

DIMFNSION    QHC(27) ,SS0RX(27) 

DIMENSION    VARXC2  7)  ,YHATHI < 27  I , YHATL0(27 ) 

INTFGFP*2    T 

D*TA     Tf/l.<J43/,T7/l.B95/,T8/l.P60/,T?/l.P33/, 
1T11/1.7C6/,T 12/1. 782 /,T10/1.B12/ 

COMPUTF  YMEAN  cor  FUTURE  USE 

DO    1    J    =1,26 
SUMY    =    CO 
K    =     I  (  J 1 
DO    2    L    *    1,K 

2  SUMY    =    SUMY    ♦    Y(JfL) 
1     YMFAN(J)     =    SUMY/IU) 

COMMUTE    BETA    FOR    Y    =    ALPA    ♦    BETA*X     FOR     EACH    HP    THE    26 
ENGINES 

Dn    3    M    =    1,26 
SUMX    =    CO 
N=!(M) 

00    A    NN    =    1,N 
4     SUMX    =    SUMX    ♦    X(M,NN) 

3  XMFAN<M)     =     SUMX     /     I(M) 
DO    5    KK    =    1,26 

SUMP    =    CO 

ssx    =  c.c 

KI     =     I(*K) 

DO    6    LL    =    1,KI 

SUMP    =    *UMR    ♦    ( (X( KK,LL)-XMEAN(KK) ) * ( Y ( KK , LL ) -YMEAN ( KK 
1)  )) 
6    SSX    =     SSX    +    (X( KK.LL)-XMEAN(KK) )**2 

S^ORX(KK)     =     SSX/U  (KK)     -    1  ) 
">    PETA(KK)    =    SUMP/    SSX 


COMPUTP    ALFA    FOP 
ENGINES 


Y    =    ALFA    +    BETA*X     EQP     EACH    0^    THE    26 


DO    12    IL=1,26 
12    ALPA(IL)    =    YMEAN(IL)     -    BET  A ( I L) *XMF AN( I L ) 
DO    7       JX    =    1,26 
XX    =    KJX) 
DO    7       KXX=1KX 
7    CftPY(JXfKXX)     =       ALCA(JX)    ♦    BETA( JX)*X( JX,KXX) 

NOW    COMPUTE    THE    ESTIMATE    OF    THE    CONOITTONAL    VARIANCE 
OP    Y    GIVEN    X    FOP    EACH    OF    THE    26    CNGTNES 

DO    P       mo    =    1,26 

NO    =     I  (MO) 

*AM    =    C.O 

DO    ~       NMO    =     1. NO 
q    caM    =    SAM    4-<Y(*0,NNO)-CAPY(MQfNNQ)  )**2 
P    *SQR(MQ)     =    SAM/(NO-2) 

roMcy-re    THr    CORRELATICN    COPEFICIFNT    BETWFEN    X    iNC    Y, 
AND    COMPUTE    THF     VARIANCE    OE    X 
SEC    P.     15P,    NAVORD    RFPORT    3369 
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DO    16 

XY    =     C 
EX    =    0 

WYr     = 


N7     = 
.C 

.c 

c 


1,26 


r 
r 
C 
r 

C 

c 
r 


FXcO«     =    C.3 
WY-SQR    =    CO 
N70    =     I(N7) 
OG     15    N77     =1,N7D 
XY    =    XY    +     X(N7,NZZ) *Y(NZ,N77) 
=  X    =    PX    +    X(N:7,N77) 
WYF    =     WYE    ♦    Y(N7,NZ7) 
EX^OR     =    EXSOR    ♦     X(N7,NZ7)**2 
15    WYFSOR    =    WYESQR    +    Y(NZ,N77)**2 

V«RX(N7)    =    (  (  I(N7)*FXSOR)-EX**2)/(  HNZJ*(  I  (N7)-l)  ) 
If    RHO(N7)    =    (N7D    *    XY    -    EX    *    WYE )  /SQPT  (  ( N  7D*EXS  QR    - 
1    EX**2)*(N7D*WYESQP    -    WYF**2)) 

WPTTF(6,100)     (NO,ALFMNO)  ,  BETA  (NO)  ,  SSQR  (  NO  )  ,  RHO  (  NO) , 
1N0=1,26) 
ICC    FOPMATf  /////30X, 'ENGINE    NO,  •  ,ieX,'ALFA»  ,  10X  ,  '  BET  A'  ,  1GX 
lf  »VAP IANCE  •  »10X,  •RHO»,////(35X,I2,SX,F10.5,5X,F1G.5, 
?5X,F10.5,5X,F10.5) ) 
CALL    BAFTLT(*SOR,I ) 

NOW    WE    FSTABLISH    A    9C    P^R    CENT    PREDICTION    INTERVAL    FOR 
Y    FOR     FACH    OF    THE    26    ENGINES,    CONSIDERING    A    HYPOTHET- 
ICAL    X-PEAOING    OF     X    =    15C 
FOR    FACH    ENGINF.     SEE    P.     163,NAV0RD    REPORT     3369 


DO    18 
PSD    = 
K  I(  IX) 
IF(  K  I 


IX    =    1,26 

SORT(  1.0    +    1.0/KIX) 

-1)*VARX( IX) ) ) 


+  ((150. C    -    XMFAN( IX)  )**2)/( 


2C 
21 
22 
23 
24 
25 

26 

3C 


K 
K 
K 
K 
If 
K 


52 


51 
1« 


!<=( 

IF( 

IF( 

TF( 

!«=( 

Ic( 

T=T6 

GO    TO 

T  =  T7 

GO    TO 

T=T8 

GO    TO 

T=T<? 

GO    TO 
T  =  T10 
GO    T0 
T=TH 
GO    tq 
T  =  T12 
YHAT    = 
YHA  th  I 
YHA  TLO 
I'M  YHA 
I  T<>  UNC 
YHA  TH  I 
GO    TO 
YHATHI 
TP(  YHA 
JTPIJN'C 
YHA  T!  0 
Ic( YHA 
GO    TH 
YHATLO 
CONTIN 
IF(NIJV 
Tc(  Nl)M 
IP(  Nl!M 

TC(MIIM 


.EO.  P) 
.EO.  °) 
.E0.1G) 
.EO.ll) 
.E0.12) 
.E0.13) 
.E0.14) 


GO 
GO 
GO 
GO 
GO 
GO 
GO 


TO 
TC 
TO 
T0 
TO 
TT 

TG 


2C 
21 
22 
23 
24 
25 
26 


3C 
3C 
30 

30 
3C 

3C 

AL 
(IX 
(TX 
THI 

(IX 
52 
(  IX 
TLO 

(IX 

TLO 

IP 

(  IX 

UF 

el 

EL 

Fl 

cl 


FA(IX)    ♦    BFTA(IX)     *    150.0 

)    =    YHAT    ♦    T*    SORT(SSORdX)  )* 

)    =    YHAT    -    t*    SORT(SSQRdX)  )* 

(IX)     .LE.    CO)    GO    TO    c0 

YHATHI (IX) 

)    =    ITRUNC    ♦    1 


RAD 

RAD 


)    =    0.0 

(IX)     .LE.    CO) 

YHA TLO (IX) 

)    =    JTRUNC    -    1 

(IX)    .LT.    O.C) 


GO    TO    51 


YHATLO( IX)    =    0.0 


)    =    O.G 


.EQ. 
.CQ. 
.EO. 
.EO. 


1) 
2) 
3) 


GO 
GO 
GC 
GO 


TC 
TO 
TC 

T0 


301 
302 
303 

304 


56 


3C1 

<>:  1 

3C2 

4C2 

3C3 

4C3 

3CA 
40 

3C5 

4C5 

3C6 

<*ce 

3C7 
i0  7 

3CP 
40  P 

2C0 


IF( 
Ic( 
IF< 
?F( 
WRI 

FOP 

an 
wr  T 
cqp 

go 

WRI 
POP 
GO 
WRI 
PHP 
GO 
WRI 
FOP 
GO 
WRT 
FDR 
GO 
WRI 
POP 
GO 
WRI 
POP 
WR  I 
POP 
1INP 
20F 
3«RR 
4FNP 
RET 
END 


NUME 

N1JMF 
NUMF 

NU*E 
TF(  h 
M£  T( 
TO  1 
T(  f 
M£T( 
TO  1 
TF(  6 
MAT< 
TO  I 
TF(6 

TO      I 

TF(6 

MAT( 

TO     1 

TE(6 

MAT( 

TO    1 

TF(6 

m«t( 

TO     1 

TF(6 

MA  T( 

TF(6 

MAT( 

IVTO 

THE 

FCIC 

•»lc 

UPN 


I  .F 
L  .F 
I  .E 
L  .F 
t*Ci 
•1', 

?<; 

,*02 

•L', 

cc 

,4C3 

'IS 
cq 

Cq 

,4C5 

•IS 

Cc 
,*C6 

•1', 
cc 

,4C7 
CC 

,4CS 
'IS 

,2CC 

UAL 

I  MOP 
Tpn « 
X,  »H 


Q. 

0. 

0. 

0. 

) 

///// 


c) 
6) 
7) 
P) 


GO  tc 
GO  T0 
GO  tc 
GO    T 


3  OS 
3C6 

3C7 

3cB 


//57X,«0AT*    POP    ALUMINUM*) 

) 

///////E^X, 'OATA     FCR     IRTN') 

) 

///////•57X,  «D4Tft    FOR    CHROMIUM1) 

) 

///////57X,'DATA    FOP    SILVER*) 

) 

///////E7X,  'DATA    FOp    MAGNFSIUMM 

) 

///////57X,«DATA  FCR  NICKCL') 

) 

///////57X,  »DATA  CCP  COPPER1) 


) 
///// 

MNP, 
///// 

V&LUF 

DENDE 
,13X, 

I  OH  E 


f/57X,»DATA  FOP  SILICON1) 

YHATLC(NR)  ,VHATHI(NR) ,MR=1, 26) 

3<?X,  'PREDICTION    INTERVAL    FOR    AN    », 

OF    Y    FHR    A • ,//33X, 'HYPOTHET ICAL    VALUE 

NT  VARIA8LF,X  =  150 • ,////63X , 

•PRFOICTED • , /4C  X, 'ENGINE    NO, • , 14X. » LOW 

N0I,///K4X,I2,17X,FP.5,13X,F9.5)) 


C 

r 

c 
c 
I- 
c 
c 
c 
c 
c 


SUBROUTINE    BARTLT(StI) 

THIS    SUBPROGRAM   COMPUTES    BARTLCTTS    STATISTIC    FCR 
TESTING    THE    HOMOGENEITY    OF    VARIANCES    OF    t    SET    OF    K 

OBSERVATIONS.    THEAROUNENTS    ARE S     IS    THE    SET    OF 

ESTIMATED    VARIANCE*    FCR    EACH    OF    THE    SETS    OF    DATA,     AND 
I     IS    THF    NUMBER    OF    OBSERVATION^    WITHIN    EACH    SET 
OF    DATA. 

DIMENSION    SI27)  ,1(27) ,DF(27) 

INTEGER*2  I 

COMPUTF  THC  DEGREES  OF  fREpDOM  FOR  EACH  cSOR 

DO  1  J  =  1,26 

1  DFl J)  =  I< J)-2 

COVUTP    THE    TOTAL    DEGREE1;    Cc    FREEDOM-    p=F.    PAGE    22% 
BROWNLEE 

F    =    r  •  C 

00    2  "k    =    1,26 

2  F    =    F     4    DF(K) 

Compute  SAMPLE  ESTIMATF  OF  SIGMA  SQUARE 

cum  =  0,0 

DO  3  L  =  1,26 

3  SUM  =  <UM  +  (DF(L)  *  S(D) 
Sir,SOR  =  sum/f 

B  =  CO 

DO  u    M=  it26 

IF(  *<  ")  .IE.    0.0)  GO  TO  ^ 


57 


r    =    r+    PF(M)    *    *LGG(?(M) /SIGSQR) 
A    CHMTINUr 

8    =     (-B) 

CL!°F     =    CO 

HO    5    N    =    1,26 
c    PLIPF    =    FLTPF    ♦     (l.C/DF(N)-l.C/F) 

C    =    1.0    +    (  l.C/75.0)*F|_IPF 

*TAT    =    B/f 

C  BAPTLETTS    STATISTIC     IS    DISTRIBUTED    APPROXIMATELY    A* 

C  CHI     SQUA&E    WITH    25    DEGREES    OF    FREEDOM 

WPITF(6,3CO)     STAT 
30f    FO°MAT( /////53X, 'BAPTLFTTS    STATISTIC • »//57X , FIO .2 ) 

WRITF(6,3C1) 
3d    PnpMAT(/////50X,»CHI     SQUARE    STATISTIC    AT    .05    LEVEL1,// 
163X,  '37.7M 
IFfSTAT    .GT.  37.7)     GO    TO    302 

WOITF( 6,303) 
303    FORMAT! ///40X, 'ACCEPT    THE    HYPOTHESIS    OF    HOMOGENEITY    OF 
1    VARIANCES') 
GO    TO     6 
3C2    WRITF(6,3CA) 

3C4    F0RMAT(///40X, 'REJECT    THE    HYPOTHESIS    OF    HOMOGENEITY    OF 
1    VARIANCES') 
6    RFTUPN 
SNP 
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OAT?.    POP    ALUMINUM 


NGINC  NT. 

&LpA 

1 

13.3u434. 

2 

6.842*5 

3 

16.97c;*5 

£. 

0.42757 

5 

15.166*4 

6 

13.33934 

7 

10.42  277 

p 

11.34943 

Q 

*. 34820 

10 

0.85517 

11 

11.59052 

12 

20.61287 

13 

12.08619 

14 

18.33127 

1* 

18.17804 

16 

15.07593 

l7 

6.42310 

18 

11.14837 

lc 

9.73746 

20 

10.34900 

21 

18.35741 

22 

6.34110 

23 

1^.30438 

2* 

6.537^3 

25 

18.72218 

26 

7.40204 

«"=TA 


J.01470 
-0.0G841 
0.0L576 
O.OC  5  2* 
0.02224 
O.Ot  746 
0. 01707 

O.Ot  76' 
0.0* 577 
0.0tl8° 

0.02414 
0.05316 
0.0c58O 
O.OC  7 '-3 
-0.0141C 
-0.03457 
0.08271 
0.01539 
0.014*9 
O.OC 272 
0.0  3504 
0.02662 
0.01323 
0.02333 
■O.OC  968 
0.01.46  5 


VARIANCE 


28.47899 

3.1.6  6  88 
16.77626 
16.31842 
30.Q0223 
^7.84  80  5 
27.71808 
14.41651 
33.21201 

0.81777 
2?. 62318 
78.7C721 
15. 074Q1 
31.08417 
48.7352C 
30.71198 

3.16762 
10.G6873 
15.26734 
12.4C99 
8*. 08536 
17.39737 
24.53^33 

1.68822 
4C. 54111 

1.6C781 


pun 


0 

•0 
C 
n 


1  =  321 
37964 
11*72 
12630 
33111 
07*15 
31Q05 
13D60 

71203 

0*885 
35753 
23280 
822  =  3 
177QQ 
07127 

SR3QA 

94023 
5C0  5  6 
0.40384 
0.0*320 
0.31535 
0.41010 
16335 
39707 
15742 
3  3606 


v; 
0 

0 

0 

0 
0 
0 
0 

0 

.  t 

•  p. 

0 
0 


BARTLETTS    STATISTIC 
76.73 


CHI     SQUARE    STATISTIC    at     .05    LEVEL 
37.7 

REJECT    THE    HYPOTHESIS     HE    HOMOGENEITY    PF    VARIANCES 
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PATA    FOR     IRON 


NGINf     NO. 

*LFA 

BETA 

VARIANCE 

OHO 

1 

5C. 20422 

0.11448 

114.27705 

0.60787 

2 

43.^6326 

-0.02378 

44.85501 

-0.29018 

3 

52.85899 

0.07885 

62.23904 

0.63434 

L 

4^.46271 

0.06691 

42.35352 

0.70755 

5 

57.84926 

0.19398 

1^7.84903 

C. 83873 

6 

65.88390 

0.03499 

693.29321 

0. 13161 

7 

63.639->6 

0.21038 

765.40845 

0.6C084 

8 

175.39531 

-0.3C968 

314.39453 

-0.75280 

Q 

44.46706 

0.191J3 

242.69214 

0.84279 

10 

46.93655 

0.05693 

65.14589 

0.49748 

11 

59.78769 

0.0530C 

86.07813 

0.380*5 

12 

3°. 87607 

0.35560 

121.76297 

C. 78992 

13 

48.83186 

0.22904 

184.28542 

0.86796 

1* 

77.23982 

-0.00702 

521.75610 

-0.04170 

1* 

41.754^3 

0.21654 

342.75439 

0.38240 

16 

63.06525 

-0.10498 

205.27632 

-0.64542 

17 

45.61081 

0.3C287 

24.24014 

0.96454 

18 

39.40042 

0.12059 

88.71730 

0.84245 

19 

57.02637 

0.06641 

129.07980 

0.56075 

20 

61.74995 

0.02155 

99.17853 

0.22816 

21 

79.81621 

0.10  270 

597.93677 

0.35013 

22 

53.77707 

0.05946 

69.42007 

0.44931 

23 

55.58455 

0.12508 

311.300  54 

0.41282 

24 

36.29274 

0.0P333 

66.22913 

0.75678 

25 

84.57494 

-0.01996 

731.59912 

-C.C8520 

26 

37.28738 

0.02948 

29.57013 

0.46688 

BARTLETTS    STATISTIC 
81.83 


CHI    SQUARE    STATISTIC    AT     .05    LEVEL 

37.7 

REJECT    THE    HYPOTHESIS    OF    HOMOGENEITY    OF   VARIANCES 
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DATA  FOR  COPPER 


NRINE  NO, 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

IS 

20 

21 

22 

23 

24 

25 

26 

#LFA 


16.10133 
7.29360 
5. 36026 
7.51995 

11.68334 
8.48832 
9.83609 
8.07563 
2.78638 
7.53345 
9.31796 
4.03100 

11.09075 

16.83313 
5.09610 

14.84856 
3.32  089 
9.26964 
8.68514 

11.42535 

12.43534 
7.85786 
5.55491 
9.97796 

16.07785 
7.91086 


BETA 

VARIANCE 

RHO 

0.01517 

28.13663 

0.20029 

0.01121 

3.33279 

0.46456 

0.06126 

2.28625 

0.95766 

0.04244 

7.06022 

0.8*119 

0.0*901 

13.07191 

0.8G406 

0.01084 

31.99890 

0.18802 

0.03832 

24.0C394 

0.61076 

0.06672 

10.72700 

0.80012 

0.08577 

24.03470 

C. 91274 

0.01489 

3.99*20 

0.51817 

0.02058 

7.38620 

0.47882 

0.06392 

2.42263 

0.85400 

0.04996 

18.08391 

0.77262 

0.0G636 

57.65799 

-0.11307 

0.05363 

7.54217 

0.56845 

•0.01592 

14.05688 

-0.43984 

0.13100 

3.68555 

0.97090 

0.02552 

4.68321 

0.82145 

0.02103 

28.38513 

0.41594 

0.01353 

11.07911 

0.40231 

0.02057 

11.64888 

0.47266 

0.01671 

15.25947 

0.28359 

0.01552 

5.26603 

0.39678 

0.03737 

9.29325 

0.81096 

0.01059 

50.29272 

0.17058 

0.02306 

5.28834 

0.69873 

BARTLETTS    STATISTIC 
72.41 


CHI    SQUARE    STATISTIC    AT     .05    LEVEL 
37.7 

PEJECT    THE    HYPOTHESIS    OF    HOMOGENEITY    OF    VARIANCES 
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OATft    FOR    <ILVER 


cisiniNF   mo. 

ALFA 

1 

3.22107 

2 

0.91387 

3 

8. 30  336 

L 

6.01645 

5 

6.18161 

6 

7.14782 

7 

5.00347 

3 

2.44647 

Q 

3.085P6 

1C 

3.12557 

11 

6.60215 

12 

2.24322 

13 

5.Q2  650 

1^ 

9.48374 

IS 

-0.02  35  5 

16 

8.42144 

17 

-0.56630 

18 

6.49206 

1° 

2.57543 

20 

7.36988 

21 

3.94133 

22 

1.12335 

23 

1.22760 

2^ 

3.38863 

25 

8.50132 

26 

2.92798 

RET  A 


0.01821 
0.0C548 
0.00955 
0.0C471 
0.02642 
0.0C980 
0.02554 
0.02110 
0.03668 
0.01790 
-0.00961 
0.02413 
0.01843 
-0.00355 
0.0*886 
-0.02204 
0.06101 
0.01983 
0.0224* 
0.00182 
0.01677 
0.0135° 
3.00936 
0.01132 
0.0C136 
-0.0C034 


VAPIANCE 

RHO 

7.22279 

0.43585 

1.42715 

0.36451 

4.05259 

0.36275 

4.70592 

C. 20703 

6.30180 

0.72415 

20.69699 

0.21049 

17.9*392 

0.51101 

6.06984 

0.48917 

12.66748 

0.79619 

2.48596 

0.67819 

1.52119 

-0.48929 

2.18181 

0.54682 

2.22144 

0.78825 

10.91075 

-0.14432 

5.32845 

0.59947 

7.28558 

-0.68556 

4.77864 

0.85632 

3.50315 

0.79120 

12.42607 

0.59432 

7.18613 

0.07334 

4.41735 

0.57907 

6.06122 

0.36266 

3.55870 

0.30231 

3.61304 

0.55839 

18.67505 

0.03652 

4.1^473 

-0.01628 

BARTLETTS    STATISTIC 
47.60 


CHI    SQUARE    STATISTIC    AT     .05    LEVEL 
37.7 

RFJECT    THE    HYPOTHESIS     OF    HOMOGENEITY    OF   VARIANCES 


62 


DATA    F0°    fHRCMIMM 


ENGINE    NO. 

ALCA 

RFTA 

VAD IANCF 

phh 

1 

-2.44742 

0.0*7  30 

25.434  13 

0.69028 

2 

-0.97^5 

0.04212 

9.3°0  5C 

0.76127 

3 

13.803B4 

-0.0*601 

17.7C932 

-0.663LC 

6 

7.7Q363 

-0.01336 

24.9  5369 

-0.25200 

c 

7.61532 

O.OC  654 

17.21991 

0.15536 

6 

14.70743 

-0.01917 

21,98987 

-0.37304 

7 

14.52510 

-0.03000 

42.90^73 

-0.41159 

a 

6.42787 

-O.CC  533 

63.19046 

-0.04386 

9 

4.43603 

0.02100 

34.7*403 

0.41414 

1C 

-0.54279 

0.06128 

18.2C578 

C. 75910 

11 

7.92  736 

0.0*316 

136.4P221 

0.25714 

12 

11.38564 

-0. 0^720 

85.0340° 

-0.27968 

13 

6.93^63 

0.03CS3 

22.8P821 

0.55148 

1* 

12.97217 

-O.Ot  3  34 

51.34  3  6° 

-0.15603 

1^ 

6.83300 

0,04071 

23.67117 

0. 2833 3 

16 

10,59030 

-0.02532 

22.68594 

-0.53017 

IT 

2.4G7Q4 

0.07504 

2P.04414 

0.6440C 

IB 

2.73783 

0.01266 

16.90269 

G. 35190 

10 

4.06411 

0.0110C 

25.95224 

o. 24264 

20 

3.93  629 

0.00669 

21.74173 

0.  15354 

21 

12.25228 

0.0  29  87 

38.51898 

0.39367 

22 

0.°2  340 

0.02031 

16.82«33 

0.33663 

23 

=.61597 

0.02869 

33.8P646 

0.300  A  9 

24 

4.71603 

0.0(325 

23.44861 

0.075  5  8 

25 

7.11337 

-O.OC  143 

29.17253 

-0.03066 

26 

7.21SQ3 

-0.01344 

1°. 43361 

-0.28456 

BARTLETTS    STATISTIC 
30.91 


CHI     SOl'APP    STATISTIC    AT     .05    LEVEL 
37.7 

ACCEPT    THE    HYPOTHESIS    OF    HOMOGENEITY    OF    VARIANCES 
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Ktjfi    FOP    MAGNESIUM 


NGINF     MC. 

ALFA 

BETA 

VAP IANCE 

PHP 

1 

1.27106 

0.01739 

0.38782 

0.64673 

2 

0.5133* 

-0.01203 

0.22000 

-0.3^639 

3 

4. 02446 

0.0C1  )7 

C. 28629 

0.28°62 

' 

2.59075 

-0.0L787 

0.39171 

-0.77457 

p 

2.27420 

9.0G84C 

C.38G63 

0.66619 

6 

2.96552 

-0.0C121 

3.11676 

-0.06816 

7 

0.58394 

0.0C5  22 

1.44880 

0.3°353 

p 

4.4515° 

-0.0C715 

1.42679 

-0.36431 

c 

1.63911 

0. Of  451 

2.23991 

0.35906 

1G 

0.58457 

-0.0C021 

0.31716 

-0.03084 

11 

1.27163 

0.0C512 

0.67668 

0.40986 

12 

1.16667 

0.0C094 

0.9*859 

O.G3950 

13 

3.86771 

0.02291 

2.310  24 

0.8^098 

14 

2.89661 

0.00057 

1.99392 

0.05512 

15 

0.29080 

0.02112 

0.52083 

0.71935 

16 

1.66836 

-0.00554 

0.76587 

-0.5899C 

17 

2.04Q64 

0.01263 

0.24931 

0.83256 

IP 

1.32456 

0.0C409 

1.04165 

0.43995 

ic 

0. 84 127 

0.0L396 

0.51687 

0.53797 

2C 

u. 50504 

0.0C237 

0.93411 

0.25669 

21 

2.41657 

0.0C083 

0.66119 

0.09060 

22 

0.46174 

0.01080 

C. 67101 

0.6  805  8 

23 

0.66353 

0.01223 

2.05501 

0.478°3 

24 

1.39193 

-0.0C162 

0.44549 

-0.26511 

25 

3.7Q430 

-0.00276 

2.01043 

-0.21981 

26 

2.22691 

-0.0C869 

0.76340 

-0.69563 

BARTLETTS    STATISTIC 
46.66 


CHI     SQUARE    STATISTIC    AT     .05    LEVEL 

37.7 

"EJECT  THE  HYOOTHPSIS  OF  HOMOGENEITY  OF  VARIANCES 


64 


DAT^     coo    NICKEL 


NT, I  MP  NO. 

ALFA 

1 

0.28900 

2 

-0.00302 

3 

1.27335 

A 

1.53143 

5 

C. 68643 

A 

1.30414 

7 

3.C0174 

8 

2.6856A 

o 

0.50666 

1C 

0.04  801 

11 

1.54.203 

12 

1.21385 

13 

1.37346 

1* 

2.34^92 

15 

-1.48613 

16 

2.87255 

17 

2.0744A 

18 

0.80567 

lc 

0.97000 

2C 

1.22195 

21 

3.03  327 

22 

1.29714 

23 

1.90  540 

2* 

C. 11*14 

25 

2.32527 

26 

1.13216 

PCTA 


0.01329 

0. Of  387 

O.0C083 

■0.00261 

0.01023 

0.0C0D7 

-0.0CQ79 

•0.00787 

0.0C704 

0.0102° 

-0.01  175 

■0.  Of.  60  2 

O.OC  83K 

-3.01  390 

0.0  23  54 

■0.0C79R 

-0. Of  434 

0.OC2O5 

0.0C015 

0.0t373 

■0.00404 

0.00122 

■0. Of  371 

0.0C766 

-0.0C057 

■O.OC  246 


VARIANCE 


4.09743 
0.55725 
1.47447 
1.74781 
0.77005 
1.73732 
2.40663 
2.06767 
2.7462* 
3.33255 
1.23409 
0.^4202 
1.14372 
1.64065 
0.73462 
1.28040 
u. 56314 
0.54576 
0.888  59 
1.63718 
1.51671 
0.65921 
2.27705 
0.89592 
7.69561 
0.85947 


RM1 


0.42437 

0. 40512 

•0.05634 

•0.18892 

0.758  38 

O.OC 5 37 

•0.04995 

•0.33742 

0.4770  8 

u.  41644 

•0.  11281 

•0.24079 

0.62830 

•0.38975 

0.69697 

•0.63039 

■0.35782 

0.42537 

0.01831 

0.30075 

•0.28009 

0.1057  2 

0.15531 

0.67501 

•0.02387 

•0.24979 


BARTLETTS    STATISTIC 
42.43 


CHI    SQUARE    STATISTIC    AT     .05    LEVEL 
37.7 

REJECT    THE    HYPOTHESIS    OF    HOMOGENEITY    OF    VARIANCES 
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DAT*     FOR    SILICON 


ENGINF    NO. 

ALFA 

BETA 

VARIANCE 

PHO 

1 

0.97502 

0.02269 

5.61894 

C. 56473 

2 

3.56379 

0.0C047 

1.31088 

U. 03513 

■* 

R.OS«?A 

0.0C02C 

1.14555 

0.01^65 

L. 

3.375*4 

-0.0(40  2 

1.43C09 

-0.3111C 

^ 

o. 544Q2 

0.00392 

^5. 180 30 

C.C4J  ;4 

6 

2.75572 

0.0C857 

3.9C288 

C. 27^89 

7 

2.99024 

0.0C626 

2.60247 

0.35762 

8 

3.12190 

0.0C500 

5.37198 

0.13986 

Q 

5.23*63 

-0.0005* 

7.83646 

-0.02473 

10 

1^.02483 

-0.05174 

54.30548 

-0.49574 

11 

2.45637 

0.OC851 

2.77044 

0.34544 

12 

6.24562 

0.00145 

18.00914 
16.C4675 

0.01363 

13 

3.48370 

0.00775 

0.19656 

14 

4.74273 

-0.00097 

1.43503 

-0.1C929 

15 

2.12191 

0.020  54 

0.90045 

C.6C832 

16 

3.54305 

-0.0C246 

2.20525 

-0.18763 

17 

1.705Q3 

0.01674 

C. 49007 

0.81773 

18 

3.12749 

0.00561 

1.89773 

0.44513 

19 

3.28432 

0.0C311 

10.41553 

0.11099 

2C 

2.01364 

0.00472 

3.13833 

0.27742 

21 

4.27623 

0.0C347 

1.21687 

0.26992 

22 

2.61480 

0.00685 

1.68334 

0.34890 

23 

2.22112 

0.01514 

3.87579 

0.44118 

24 

2.34076 

0.00553 

4.85940 

0.27276 

25 

4.16613 

0.00015 

3.94967 

0.00901 

26 

3.82401 

-0.00690 

3.46491 

-0.33962 

BARTLETTS    STATISTIC 
158.51 


CHI    SQUARE    STATISTIC    AT    .05    LEVEL 
37.7 

REJECT    THE   HYPOTHESIS    OF    HOMOGENEITY    OF   VARIANCES 
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APPENDIX  H 
GENERAL  RESULTS  FOR  JET  ENGINES 

Essentially  the  same  analysis  was  conducted  for  jet  models 
J52-P-6,  J60-P-6,  and  J33-A-2M-,  as  was  conducted  for  model  R2000-6. 
When  scatter  diagrams  for  the  jet  models  were  plotted  they  pro- 
vided no  encouragement  that  a  functional  relationship  might  exist 
between  the  independent  and  dependent  variables.   As  part  of  the 
analytical  analysis  a  stepwise  regression  was  executed  in  order  to 
determine  which  of  the  two  variables,  hours  since  engine  oil  change 
or  hours  since  engine  overhaul,  was  the  better  predictor.   For  all 
three  models  the  results  were  inconclusive .   Approximately  one- 
half  the  engines  showed  hours  since  engine  oil  change  to  be  the 
better  predictor  and  for  the  other  half  hours  since  engine  overhaul 
was  the  better.   Simple  regressions  were  executed  using  each  of 
the  two  independent  variables.   Regardless  of  which  variable  was 
used  in  the  equation,  about  half  of  the  regression  coefficients 
were  negative.   This  meant  that  in  about  half  of  the  jet  engines 
there  was  a  decrease  in  the  estimate  of  the  density  of  wear 
metals  as  time  since  oil  change  increased.   Inconsistent  behavior 
such  as  this  could  only  be  interpreted  as  evidence  that  neither 
of  the  two  control  variables  possessed  reliable  predictive  powers 
for  the  dependent  variable. 

For  each  model  Bartlett's  test  for  homogeneity  of  variances 
was  conducted.   The  results  were  similar  to  those  obtained  for 
the  reciprocating  engine;  i.e.,  the  jet  engines  of  a  particular 
model  could  not  be  considered  as  coming  from  the  same  population. 
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Therefore,  as  with  the  reciprocating  engines,  data  manipulation 
must  continue  to  be  carried  out  separately  for  each  engine. 

When  the  raw  data  for  jet  engines  was  compared  against  raw 
data  for  a  reciprocating  engine  an  immediate  difference  was  detec- 
ted in  the  densities  of  all  elements.   The  densities  in  parts  per 
million  for  jet  engines  was,  in  general,  considerably  less  than 
that  for  a  reciprocating  engine.   A  typical  comparison  is  shown  in 
Table  II.   For  model  J33-A-24  five  engines  were  analyzed,  for 
model  J52-P-6  twenty-four  engines  were  analyzed,  and  for  model 
J60-P-6  twenty-one  engines  were  analyzed.  No  engine  was  included 
in  the  analysis  unless  it  had  eight  or  more  data  points.   Since 
the  estimated  densities  for  jet  models  are  of  the  magnitude  shown, 
it  is  entirely  possible  that  these  densities  are  within  the  noise 
level  of  the  spectrometer.   This  probably  explains  the  fact  that 
neither  hours  since  engine  oil  change  nor  hours  since  engine  over- 
haul appear  to  be  of  value  in  predicting  the  internal  condition  of 
engines  chosen  from  the  set  of  three  jet  models  analyzed  in  this 
study. 
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Density  of  Element  (PPM)  by  Model 

Element 

R2000-6 

J60-P-6 

J52-P-6 

J33-A-24 

Aluminum 

7-25 

0-2 

0-5 

0-4 

Iron 

58-100 

1-5 

3-6 

1-3 

Chromium 

3-18 

1-7 

1-9 

1-8 

Silver 

5-11 

0-1 

0-2 

0-3 

Copper 

9-26 

0-1 

0-1 

0-2 

Magnesium 

2-4 

0-4 

0-4 

0-1 

Nickel 

0-6         0-2 

1 

0-2 

0-7 

TAI 

ILE  II 

Comparison  of  Typical  Densities  Between  Jet 
And  Reciprocating  Engines 
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